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This paper deals with a cosmological model in which the universe is filled with tachyon dark energy 
in order to describe current and future accelerating expansion. We obtain that the simplest condition 
for the regime of phantom energy to occur in this scenario is that the scalar field be Wick rotated 
to imaginary values which correspond to an axionic field classically. By introducing analytical 
expressions for the scale factor or the Hubble parameter that satisfy all constraint equations of the 
used models we show that such models describe universes which may develop a big rip singularity 
in the finite future. It is argued that, contrary to a recent claim, the entropy for a universe filled 
with dark energy is definite positive even on the phantom regime where the universe would instead 
acquire a negative temperature. It is also seen that, whichever the fate of the tachyonic accelerating 
universe, it will be stable to any fluctuations of the scalar field, and that since the considered models 
have all an imaginary sound speed, any overdense regions will undergo an accelerate collapse leading 
rapidly to formation of giant black holes. Finally the conjecture is advanced that these black holes 
may be the supermassive black holes that most galaxies harbor at their center. 

PACS numbers: 98.80.-k , 98.62.Js 



I. INTRODUCTION 

The discovery that high redshift supernovae of type 
la are less bright than expected [1,2] has opened a re- 
ally new cosmological scenario with far-reaching implica- 
tions which may even reach the standard model of par- 
ticle physics and the very nature of general relativity it- 
self. The straightforward interpretation of this result is 
that the universe is currently undergoing a period of ac- 
celerated expansion, similar to primordial inflation, but 
with a rather uncertain future. A most conservative ap- 
proach to explain the result while keeping general rela- 
tivity essentially untouched is invoking the inclusion of a 
now dominating extra component with negative pressure, 
usually known as dark energy [3] , in the universe. There 
are several candidates for dark energy whose properties 
are being probed by cosmological observations, includ- 
ing mainly a pure cosmological constant, quintessence [4] 
or k-essence [5] scalar fields (which may [6] or may not 
be tracked [7]), and some unified dark matter scenar- 
ios among which two main plausible models have been 
suggested: the generahzed Chaplygin gas [8], which is 
described by a highly non conventional equation of state, 
and the cosmic tachyon theory based on non canonical 
kinetic energy [9]. 

On the other hand, the increasing dropping of de- 
tailed analysis leading to quite an ample observation- 
ally acceptable parameter space beyond the cosmological- 
constant barrier [10] is opening the really intriguing pos- 
sibility that the universe is currently dominated by what 
is dubbed as phantom energy [11]. Phantom energy has 
rather weird properties which include [12]: an energy 
density increasing with time, naively unphysical super- 
luminal speed of sound, violation of the dominant en- 
ergy condition which would eventually allow existence of 
inflating wormholes and ringholes [13], and ultimately 



emergence of a doomsday singularity in the finite fu- 
ture which is known as the big rip [14]. The regime for 
phantom energy takes place for state equation parame- 
ters oj = p/ p < —1 and has been shown to occur in all 
current dark-energy models. However, whereas the big 
rip singularity is allowed to happen in quintessence [14] 
and k-essence [15] models, it is no longer present in mod- 
els based on generalized Chaplygin-gas equations of state 
[8,16] having the form P = — A/ p^^ with A and n being 
constants. No discussion has been so far made never- 
theless on the occurrence of phantom energy and big rip 
in the other major contender model for unified dark en- 
ergy: the tachyon matter scenario of Padmanabhan et al. 
and others [9,17] or its sub-quantum generahzation [18]. 
Abramo and Finelli have in fact used [19] a Born-Infeld 
Lagrangian with a power-law potential and recovered a 
nice dark-energy behaviour, but did not consider the neg- 
ative kinetic terms which appear to characterize phantom 
energy. 

While defining a phantom energy regime in such sce- 
narios appears to be rather straightforward, it is quite 
more difficult to obtain an associated expression for the 
scale factor of the accelerating universe which allows us to 
see whether or not a big rip singularity may occur. Under 
the assumption of a constant parameter for the equation 
of state, we derive in this paper a rather general solution 
for the scale factor of a universe dominated by tachyon 
matter, and show that, quite similarly to as it happens in 
current quintessence models [2,3,6,7], that solution may 
predict both an ever expanding accelerating behaviour 
and the occurrence of a big rip singularity, theoretically 
covering therefore all observationally allowed possibilities 
other than a future big crunch. 

Unified dark matter or energy models are very inter- 
esting [8,9] in that they are able to simultaneously de- 
scribe the properties of both dark matter and some form 
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of dark energy (often a cosniological constant) as two 
opposite limiting situations (at high and low densities 
respectively) arising from just one theory. The general- 
ized cosmic Chaplygin models have been however claimed 
to be observationally inconsistent and actually unstable 
to short-scale fluctuations [20]. This conclusion was ex- 
tended to also encompass any unified dark matter or en- 
ergy models, including the cosmic tachyon scenario [20]. 
However, that conclusion is true only when both baryons 
and dark matter are assumed to vanish in the Chaplygin 
model. If a given, finite fraction of either is introduced 
in this model then the oscillations may turn out to be- 
come very small and do not necessarily spoil the spectra 
[21,22]. 

In the present paper we shall use our general solution 
for the scale factor in a universe filled with tachyon dark 
energy in order to analyze in some detail the nature and 
evolution of the fiuctuations. We obtain that the uni- 
fied tachyon theory is stable to any fluctuations and that 
these may contribute the structure in the universe in a 
rather decisive way. Actually we conjecture that the nec- 
essary accelerate collapse of the formed overdense patches 
may lead to the formation of the kind of giant black holes 
which are now believed to occur at the center of most 
galaxies [23]. 

The paper can be outlined as follows. In Sec. II 
we discuss a rather general solution for the scale fac- 
tor which satisfies all requirements and constraints im- 
posed by tachyon theory. The phantom regime for such 
a solution is then investigated in this section where it is 
also seen that it shares all funny properties of quintessen- 
tial phantom energy, including the new characteristic of 
having a negative temperature, even though its entropy 
keeps being definite positive for any value of the equa- 
tion of state parameter. The evolution of fluctuations 
with small wavelengths is considered in some detail in 
Sec. Ill, using both a synchronous perturbation formal- 
ism and the relativistic analog of the Newtonian first or- 
der perturbation equation. Sc;c. IV contains a discussion 
on the gravitational spherical collapse of formed over- 
density regions which is further accelerated by pressure 
negativencss and the very nature of the tachyonic stuff. 
The conjecture that the giant black holes formed as a re- 
sult of that collapse could actually be those black holes 
which most galaxies harbor at their central bulge is ad- 
vanced in Sec. V. Results are summarized and briefiy 
discussed in Sec. VI. 



II. TACHYON MODEL FOR DARK ENERGY 



ter, and P is the corresponding pressure. We shall re- 
strict ourselves to consider a description of the current 
cosmic situation where it is assumed that the tachyon 
component largely dominates and therefore we shall dis- 
regard in what follows the non-relativistic and relativistic 
components of the matter density and pressure. For the 
tachyon field <j) we have [17] 
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(2.2) 



in which V{^) is the tachyon potential energy. Assuming 
an equation of state P = u)p for the tachyon matter, we 
then deduce that 



w = 0^ - 1. 
Finally, the equation of motion for (j) is 

1 dV{(l)) 
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V{(j)) dcp 



(2.3) 



(2.4) 



We shall show next that there exists a general solu- 
tion for the scale factor a{t) in this tachyon-field sce- 
nario which has exactly the same dependence on time as 
that has been obtained in the general solution for a pure 
quintessence scalar field, and hence wc also show that 
such a scenario admit the existence of a tachyon phan- 
tom field which leads to a singularity in finite time. In 
fact, a recipe has been provided by Padmanabhan him- 
self [17] according to which, given the explicit form for 
the scale factor a{t), a complete specification of the full 
(/)-field theory can be achieved by using the following re- 
lations: 
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(2.6) 



(2.7) 



For the most favored cosmological spatially flat sce- 
nario, the Priedmann equations read 
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(2.1) 



where p ~ Pnr + Pr + Pep is the energy density for, re- 
spectively, non-relativistic, relativistic and tachyon mat- 



Our task then is to choose a general expression for a(t) 
which simultaneously satisfles relations (2.5), (2.6) and 
(2.7), together with the Friedmann equations (2.1) and 
the equation of motion for the tachyon field (2.4), which 
be able to match the accelerating behaviour of the current 
universe and implies a physically reasonable and suitably 
motivated field potential. If we assume a linear time- 
dependence of the tachyon field 4> and hence constancy 
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of parameter w, then it is not difficult to check that a 
general form of such an expression for a{t) can be written 
as 



a{t) = 



2/[3(l+a-)] 



(2.8) 



where ao is the initial value of the scale factor at the 
onset of tachyon dark energy domination. We note that 
this solution describes an accelerating universe in the in- 
terval — 1/3 > uj > —1. At the extreme point uj = —1/3, 
a{t) describes a universe whose size increases just as t, 
such as it should be expected. It is worth realizing that 
by simply trivially re-scaling the time parameter, solu- 
tion (2.8) turns out to be nothing but the scale factor 
that represents the most general solution for the case 
of a quintessence scalar field for a constant equation of 
state [13]. On the other hand, for a scale factor (2.8) the 
tachyon field and potential are given by 



= ^0 + VT+w t 



(2.9) 



B. w < -1 

The phantom energy regime will be characterized by 

values of the state equation parameter such that uj < 
— 1 and a consequent violation of the dominant energy 
condition, i.e. 



P + p 



T/(0)</.2 



1 



< 0. 



(2.12) 



Such a regime can be obtained by simply Wick rotating 
the tachyon field so that (f) ii>, with which the field 
4> can be viewed as an axion tachyon field [25], as the 
scale factor a{t) and the field potential V{^) keep being 
positive and given respectively by 



a{t) 



^0 



-3(|a;|-l)/2 



2lH 



l)t 



-2/[3(|a,|-l)] 



(2.13) 



which accounts for a big rip singularity at finite future 
time 
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(2.10) 



We note that as i/i ^ oo this potential reasonably van- 
ishes after taking the form already considered by Pad- 
manabhan and others [17]. As </> — > <;io at * = 0, V{(p) 
tends to a finite constant value, so clearly separating 
from the unphysical behaviour of the potential consid- 
ered by Padmanabhan and compatible with what can be 
supported by string theories [17]. We regard therefore 
potential (2.10) as being physically reasonable. Finally, 
using Eqs. (2.2), (2.3) and (2.8)-(2.9), we obtain for the 
speed of sound 



P 
P 



(2.11) 



This result is in contradiction with the value c\ = —to 
considered in Refs. [20,21]. The reason for that discrep- 
ancy comes from the definitions used in this paper of 



the Lagrangian [17], L = P = —V{4>)^J\ — (j)^, and the 
energy density given in Eq. (2.2). 

For the accelerating-expansion regime, we see thus that 
the speed of sound becomes imaginary, a case which could 
imply an accelerating collapse of the tachyon stuff that 
can still be circumvented however [24]. Eq. (2.11) is 
compatible with all the above requirements, provided we 
assume a linear time-dependence for the field 0. On the 
other hand, since w < 0, this equation may be also com- 
patible with a wave equation for tachyon fluctuations 
where there is factor —uj in front of the Lagrangian, such 
as the one used by Frolov, Kofman and Starobinsky [20] . 



and 

y($) 



3(|a;|-l)a^(l"l-^)/^' 



(2.14) 



87rG 



-3(|u;|-l)/2 3 



2 ' 



.0 - - ^o) 

(2.15) 

with $0 — «^o- We note that both this potential and 
the phantom tachyon energy density. 
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(which has been obtained by using Eqs. (2.2), (2.3) and 
(2.15)) increase with time up to blowing up at the singu- 
larity at t = t*, to steadily decrease toward zero there- 
after, in a regime where the universe contacts to finally 
vanish as f ^ oo. Thus, the tachyon model for dark 
energy contains a regime for phantom energy which pre- 
serves all the weird properties shown by this in current 
quintessence and k-essence scenarios; i.e. superluminal 
imaginary speed of sound, increasing energy density, vi- 
olation of dominant energy condition and a big rip sin- 
gularity followed by a contracting phase. 

It has been quite recently pointed out [26] that the en- 
tropy for phantom energy is definite negative, and that 
therefore a cosmic regime with phantom energy should 
be excluded. However, there exists a rather general ar- 
gument, valid for all phantom models studied so far in- 
cluding the one considered in this paper, which appears 
to prevent these conclusions. In fact, for a general FRW 
flat universe filled with a dark energy satisfying the equa- 
tion of state p = Lup (with cij=const.) where the first law 
of thermodynamics holds and the entropy per comoving 
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volume stays constant [27], the temperature of the uni- 
verse is given by 



T = k(1 + w)a 



-3a; 



(2.17) 



where k is a positive constant. From Eq. (2.17) it is 
aheady seen that T is negative for w < —1. On the other 
hand, it is well known that by integrating the cosmic law 
of energy conservation, p + 3p(l + lo)H = 0, we have 



Pod 



-3(l+a>) 



(2.18) 



Now, from Eqs. (2.17) and (2.18) one can readily derive 
the following generalized Stefan-Boltzmann law 



k(1 + w) 



(2.19) 



Note that for uj = 1/3 Eq. (2.19) consistently reduces 
to the usual law for radiation, and that for > a; > —1 
the dark energy density decreases with the temperature, 
which is already a rather weird behaviour. However, for 
the regime where uj < —1, in order to preserve p positive, 
we must necessarily take T < 0, and therefore p will in- 
crease with |T|. Since on the phantom regime p increases 
with the scale factor a{t) it also follows that |T| increases 
as the universe expands on that regime. Finally, a gen- 
eral expression for the dark energy entropy can also be 
obtained [26] which reads 



S = Co 



T 



1 



i/c 



V, 



(2.20) 



where Cq is a positive constant and V is the volume of 
the considered portion within the dark energy fluid. It 
follows that, contrary to the claim in Ref. [26], the en- 
tropy of a dark-energy universe is always positive, even 
on the phantom regime. Actually, by inserting Eq. (2.17) 
into Eq. (2.20) one attains that S'=Const. when we take 
V as the volume of the entire universe. 

It is nonetheless the temperature which becomes nega- 
tive instead of entropy for a; < —1. Even though it is not 
very common in physics and therefore can be listed as just 
another more weird property of the phantom scenario, a 
negative temperature is not unphysical or meaningless. 
Systems with negative temperature have already been 
observed in the laboratory and interpreted theoretically. 
In the case of phantom energy it means that the entropy 
of a phantom universe monotonically decreased if one 
would be able to add energy to that universe. Hence, 
a a; < — 1 universe would always be "hotter" that any 
ui > —1 universe, and if two copies of the universe were 
taken, one with positive and other with negative temper- 
ature, and put them in thermal contact, then heat would 
always flow from the negative-temperature universe into 
the positive-energy universe. It could yet be argued that 
negative temperature is a quantum-statistical mechanics 
phenomenon and therefore cannot be invoked in the clas- 
sical realm. However, a negative temperature given by 



Eq. (2.17) when a; < — 1 [27] can still be heuristically 
interpreted along a way analogous to how e.g. black hole 
temperature can be interpreted (and derived) without us- 
ing any quantum-statistical mechanics arguments; that 
is by simply Wick rotating time, t — » it, and checking 
that in the resulting Euclidean framework t is periodic 
with a period which precisely is the inverse of the Hawk- 
ing temperature [28] . Thus, the Euclideanized black hole 
turns out to be somehow " quantized" . Similarly in the 
present case, the phantom regime can be obtained by 
simply Wick rotating the classical scalar field, ^ — > 
which, by Eq. (2.9), is equivalent to rotating time so that 
t IT, too. It is in this sense that the phantom fields 
are also "quantized" and that the emergence of a negative 
temperature in the phantom regime becomes consistent. 



III. RENAISSANCE OF THE UNIFIED 
TACHYON DARK MATTER MODEL 



The biggest problem which has been claimed to be con- 
fronted by tachyon theory of dark energy refers to its 
property of being a unified dark matter model with a 
nonzero speed of sound. In fact such a kind of model has 
been ruled out as causing violently unphysical blowup in 
the matter power spectrum [20]. Actually, the precise 
way in which the unified dark matter models separates 
from what is observed has been only carried out for the 
case of the so-called cosmic generalized Chaplygin gas 
[20] whenever baryons and dark matter are assumed to be 
absent. Rather general arguments were then raised [20] 
to generalize the conclusion to any unified dark energy 
model whose equation of state allows for nonzero values 
of the sound speed. In what follows we shall show that, 
even though baryons and dark matter are not explicitly 
introduced, such a generalization is no longer valid, at 
least for unified models of tachyon dark matter and en- 
ergy. We will in fact see that perturbations with wave- 
length below the Jeans scale do not increase along the 
available time interval for the imaginary values taken by 
the sound speed for the tachyon field for any uj and a 
scale factor as given by Eq. (2.8). The stress tensor for 
the tachyon scalar field can be written in the perfect fluid 
form [17] 



Ti = {p + p)u'uk-p5l 



(3.1) 



where the pressure and energy density are given in Eq. 
(2.2) and 



dk<t> 



(3.2) 



This stress tensor can be now split into a pressureless, 
dark matter component and a dark energy component 

[17]: p = Pv + Pdm and p = Pv + pdm- Now, the vac- 
uum component can correspond to either a cosmological 
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constant if we choose [17] 



, Pdm = 



Pdm = 



= V{(f))\Jl - (f)^, Py = 



(3.3) 



(3.4) 



(for which cjy = —1), or to a sfowly- varying quintessential 
scalar field when we choose 



the conformal time rj, not in the physical time t, with rj 
and t being related to each other by means of Eq. (3.11). 
It can be readily seen by using Eq. (2.8) that a growing 
mode in r] necessarily implies a decreasing mode in t for 
any w < —1/3 characterizing an accelerating universe. 
We note, moreover, that the function p, is pure imaginary 
for Cg real and real for pure imaginary Cg. In the present 
model, we have from Eq. (2.11) that = w. Thus, if 
— Lo is assumed to be constant, all the problem reduces 
to study a simple differential equation given by 



PUM = V{(j))yl - 02, pdm = 



(3.5) 



2(1 -3^ 
(l + 3cf)2ry2 



At = 0. 



(3.12) 



Pv 



If the range of equation of state parameter would be 

extended to include the regime < uj < 1/3, then in 
that regime all would happen like in the usual case. 



In this latter case we have 
w„ = 1 



= ) Pv = —y{4')\l ^ ~ 4''^ ■ (3-6) that is the perturbations would oscillate for fc^c^ 



(3.7) 



which for sufficiently slowly- varying fields can even take 
on values < — 1. 



A. Fluctuation formalism 

Using the synchronous perturbation formalism first de- 
veloped by Lifshitz and Khalatnikov, the study of scalar 
perturbations, /i, around the FRW metric can be reduced 
to investigate the solutions of the equation for the fcth 
mode [29] 



p" + k'c] 



,2 2 



M = 0, 



(3.^ 



where ' = d/drj, with f] = J dt/a being the conformal 
time, cl= const, and the function p is defined as 

{h' + ajh)a i{Hh-Hh)a 
M = p = 7= > (3.9) 



with ■ = d/dt, a — a' /a = a, 7 = 1 — oi! jo? = 1 — a/d^, 
the scale factor a{i) being given by Eq. (2.8) and 



air]) = 



(1-F3a;)r? 



2/(l+3w) 



where 



1 + 3tJ 



a(t) 



(l+3w)/2 



(3.10) 



(3.11) 



At first sight it seems that the small-scale fluctuations 
must grow when < 0. In fact, taking the limit for large 
k in Eq. (3.8) we obtain that the solution should be a 
growing mode for < 0. But this is a growing mode in 



» 

2(1 — 3ti;)/[(l + 3a>)2ry2] and had an exponential behaviour 
for k^cl << 2(1 - 3tj)/[(l + 3w)2ry2]. However, if we ex- 
tend the range of cj- values to also cover values lo > 1/3, 
then gravity would behave like a repulsive force and 
played on the same team as pressure does. For the most 
interesting range where the universe shows an accelerat- 
ing expansion, lo < — 1/3, however, it is the pressure what 
plays on the same team as gravity does. In that case, if 
>> 2(l-3w)/[(l + 3w)V], then the solution to the 
differential perturbation equation reads 



p,(x e 



(/=|c?r/^.)=exp(±fc|c^|i/2 



1 + 3^ 
2^3(1 + 1..) 

1 + 3lo 



(3.13) 



with 



a, 



3(l-l-a;)/2 




+ -{l+u;)t. 



In the range of state equation parameters —1/3 > lu = 
> — 1 it can be seen that p decreases with time t, from 
an initial value 



p = po (x exp 



±2fclc2|V2 (i+3,2)/2- 

l + 3c2 "° 



(3.14) 



at t = 0, to become a minimum p (x 1, as t ^ oo. Also 
for the case where to = < —1 p starts with a value 
given by Eq. (3.14) and becomes oc 1 as one approaches 
the big rip time . That decrease of density fluctuations 
with small wavelengths clearly separates from what was 
previously predicted by Carturan and Finelli, and Sand- 
vik et al. [20], being also in sharp contrast with what 
arises from the equation that describes the evolution of 
small perturbations in Newtonian physics, that is [30], 



ei -I- [k'^cj - 47rGeo) ei = 0, 



(3.15) 



where = Pi/ei, with pi and ei being the pressure 
and energy density first order perturbations. In fact, for 
< 0, it is obtained that ei increases exponentially with 
time so that 



ei oc exp I ±W fc^jc^l -|- AnGeot 



(3.16) 
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Although the mathematical reason for having obtained 
so a different behaviour in our estimate for the time evo- 
lution of density perturbations must be addressed to the 
specific form with which the scale factor used in this 
paper depends on u (which generalizes more restricted 
possible expressions), the physical reason why the short- 
scale fluctuations are stable in our model, even in the case 
where no contribution from baryons and dark matter is 
explicitly considered, can be ultimately interpreted to be 
the unified character of that model where dark matter is 
also implicitly present. 

We shall turn next to derive general analytical solu- 
tions for the density contrast, (5, in our tachyon model, 
thereby checking and refining the conclusions attained so 
far in the above heuristic model. That study will allow us, 
moreover, to deal in a proper way with the important is- 
sue of the initial conditions for fluctuations. If we assume 
that the tachyon field is given by Eq. (2.9), the tachyon 
field theory considered in Sec. II is fully equivalent to the 
simple theory of a perfect fluid with constant equation of 
state P = uip and squared sound speed = w in a flat 
universe with a scale factor given by either Eq. (2.8) if 
ti) > — 1 or Eq. (2.13) if < —1. Thus, even though the 
density fluctuations of a perfect fluid generally behaves 
differently from the density fluctuations of a scalar field, 
in an adiabatic approach we can use the relativistic ana- 
log of the Newtonian first order perturbation equation in 
order to describe fiuctuations of the density contrast. In 
our present case, from the form of the energy density and 



the scale factor one can derive an expression for the rela- 
tivistic analog of the Newtonian first order perturbation 
equation in Fourier space [31] which reads for a density 
contrast perturbation 5k with wave vector k 



9cl)S',- £(l + 2c^-3c^)-(^^j 4 = 0, 

(3.17) 

where now ' = d/dlna and we have used ui = c^. Due to 
the presence of a first-order derivative term in Eq. (3.17), 
also in this case the taking of the simple large k limit in 
this equation when < cannot by itself uncover the 
real behaviour of small-scale fluctuations with physical 
time t. An explicit integration of this differential equa- 
tion should also be done. For this to be accomplished, 
one can now re-write Eq. (3.17) in the form 



3 



{l + 2cl-3ct)-k'cla 



dk = 0,(3.18) 



where now ' = d/da and we have used Eq. (2.8). Since 
Cg < this differential equation admits a general analyt- 
ical solution which can be expressed in terms of modified 
Bessel functions [32]. Depending on the initial condi- 
tions that solution can be taken to be either some of the 
following ones or some of their combinations 



i7r(5+3cg) 

2(l + 3c^) 



I 1 + 3c^ . 



(3.19) 



± 



2A:|. 



2(l+3cJ) V -I- + <JCs 



(l+3c^)/2 



(3.20) 



for -1/3 > c2 > -5/3, 



9 , i^(5+3cg) / 9J,| 211/2 



ITT 



2(l+3c|) 



l + 3c2 



(3.21) 



for < -5/3. 



We must now set the initial conditions for these pertur- 
bations. If we assume that dark energy is present in the 
universe from its very origin at a = 0, the simplest and 
most natural condition that such perturbations should 
satisfy is that they ought to tend to vanishing values as 
one approaches the initial singularity both for lo > —1 
and (jj < —1. One should then set Sk = as a ^ 0. In 
the limit a ^ (i.e. for arbitrarily large values of the 
argument z of the modified Bessel function, B{±z)), we 



have for B{±z) = I{±z) 

Sk oc a-(i-3-^)/2 exp ( ±iMl a(i+3e^)/2') (3.22) 
V 1 + 3c2 J 

which vanishes for the -I- sign (i.e. for — tt < argz < 7r/2) 
and blows up for the - sign (i.e. for 7r/2 < arg^ < tt), 
and for B{±z) = K{±z), 

Sk a a-(^-3^^)/2 exp UlML^a^'+^cl)^ ^ (3 23) 
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which vanishes for the lower sign and blows up for the 
upper sign. Then the above initial conditions imply that 
we have to choose the modified Bcssel functions I{+z) 
and K{—z), or some combination of them, for universes 
which start being dominated by dark energy at a = ao- 
Note that at a = 0,0, i.e. at the time when dark en- 
ergy starts to dominate, Sk consistently takes on finite, 
nonzero values. It follows that the universe is stable un- 
der perturbations 6k- Moreover, if the initial value of the 
perturbation at a = ao is small enough then the fluctu- 
ation will first grow up to a given maximum value to be 
steadily damped thereafter, vanishing as a ^ 00, both for 
Lo > —1 and < — 1. There would then be an interval of 
scale factor values where these fluctuations had observ- 
able effects. Since the sound speed is pure imaginary in 
all the considered cases, fluctuations will gravitationally 
collapse in an accelerated fashion to form, together with 
observable matter, some structures with cosmological in- 
terest. 

We note finally that consistency of using Eq. (3.17) 
with the above initial conditions is manifested, if not 
guaranteed, by the feature that they lead to exactly the 
same conclusions as those obtained from the synchronous 
model. A caveat is worth mentioning however. Once we 
have discussed fluctuations in the linear phase, it appears 
most convenient to consider the fluctuations in the non- 
linear phase, as these could become eventually dominant 
in some; rcginicis. Such a study is left for a future work. 

On the other hand, another caveat should also be 
mentioned. Even though our solutions apply to scales 
in the sub-horizon regime, if the universe is expanding 
in an accelerating fashion then it would also be natu- 
ral to consider that sooner or later all scales will be- 
come larger than the horizon so that the Newtonian ap- 
proximation would break down, specially when the limit 
t — > 00 is approached for w > — 1. In the case that 
k'^cl « 2(1 - 3cj)/[(l + auifr]^], the solutions to, for 
example, Eq. (3.12) can be approximated to 

with A and B arbitrary constants. We then see that for 
Cg < —1/3, n will grow with time t. 



B. Fluctuations in the contracting phase 

Once we have checked that our tachyon model for an 
accelerating universe is stable against density fluctua- 
tions both when w > — 1 and for w < —1 up to the big 
rip time , we next proceed to briefly consider what hap- 
pens with such fluctuations when a; < —1 at times t > t*. 
This study is of some interest for the following reason. 
Even though the big rip corresponds to a curvature sin- 
gularity that does not strictly allow for causal connec- 
tions between the regions before and after the big rip, 
this might still be circumvented by connections between 
these regions rendered physically plausible by means of 



wormholcs inflated by the accelerated expansion [13] or 
phantom energy accretion [33]. It can be seen however 
that if wo restrict ourselves to the case w < — 1 in Eq. 
(3.13) then for t > there will be a huge growth of the 
density of perturbations with small wavelengths, so ren- 
dering the evolution of the universe after the big rip fully 
unstable. On the other hand, in order for studying the 
evolution of the fluctuations after t^, in some more detail, 
one must introduce the condition that Sk ~ as t ^ 
(i.e. as a ^ 00). Now, for a —^ 00 (i.e. as z tends to 
zero), we have 

4[/(±z)] (X 0^+3'=' ^ 0, (3.24) 

Sk[K{±z)] (X 0, (3.25) 

if -1 > cl > -5/3, and 

Sk[K{±z)] cx a-3(i-'=')/2 ^ 0, (3.26) 

for < —5/3. We have to choose then the modified 
Bessel functions I{—z) and K{—z), or some given com- 
bination of them, for universes starting to evolve just 
after the big rip singularity. It can be then checked that 
such conditions also predict that the contracting universe 
that starts evolving just after the big rip is unstable to 
the considered fluctuations. 



IV. STRUCTURE FORMATION IN A 
TACHYON DARK-ENERGY UNIVERSE 

In this section we will briefly consider the fate of the 
fluctuations 6k once they have formed and evolved to 
their nearly largest size. From the onset of the epoch 
where dark energy starts dominating the formed over- 
dense regions in the universe will follow an evolution 
which is governed by a stuff characterized by negative 
pressure and speed of sound squared. That will inex- 
orably makes the force involved in such fluctuations to 
be attractive as pressure will now contribute on the same 
attractive pattern as gravity. Thus, the overdense re- 
gions would immediately start to undergo an accelerated 
gravitational collapse which involved both dark matter 
and dark energy, and hence any now subdominant possi- 
ble visible matter that could be more or less associated 
with them, towards a central singularity with infinite 
density. In practice, since the stuff making up tachyon 
dark matter-energy is expected to be spinless and charge- 
less, in the present case only some macroscopic dissipa- 
tive processes, other than pressure, spin degeneracy or 
electric charge repulsion, could intervene to convert ki- 
netic energy of collapse into random motions before the 
singularity. Black holes mainly made from dark mat- 
ter could be thereby most easily formed at the end of 
the collapse, though the residual dissipative forces could 
still stop the collapse before reaching the final dynamical 
equilibrium. 
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A. Accelerated spherical collapse 

We shall assume in what follows a spherical collapse 
model with no shell-crossing, so that wc can ignore 
the spatial dependence of the involved fields [34]. In 
principle, the evolution of a spherical overdense patch 
with scale radius R{t) containing independent dark mat- 
ter (with density pdm) and dark energy (with density 
p = p/w) is governed by the Raychaudhuri equation [30] 



3R = -AnGR [pdm + (1 + 3a;)p] . 



(4.1) 



For our unified model, the first term in squared brackets 
of this equation becomes redundant because dark matter 
is already contained in the unified tachyon dark descrip- 
tion encapsulated in the term containing density p. Thus, 
for the case being considered, Eq. (4.1) reduces to: 



3R = -4nGR{l + 3oj)p. 



(4.2) 



Now, in order to check consistency of our unified dark 
energy model, lot us first apply Eq. (4.1) to an energy 
density as given in Eq. (2.2) which, when expressed in 

terms of variable T = a^^^'^'^^^'^ + 3(1 + uj)t/2, reads 



P = 



(4.3) 



That is, one is considering the universe as a whole as a 
perturbation and, therefore, one should recover the scale 
factor (2.8) as a solution from Eq. (4.1). This is in fact 
the case, but we obtain also another solution given by 



R = ac = t(i+3'-)/[3(i+'^)1 = a(t)(i+3' 



;w)/2 



(4.4) 



Note that this additional solution will vanish as a — > oo 
both for a; > — 1 and for w < — 1 and represents therefore 
a collapsing universe. 

For a spherical overdense patch with scale radius R{t) 
within the universe, the Raychaudhuri equation will be 



3R = -4ttGR{1 + 3w){p + p), 



(4.5) 



where p is given by Eq. (4.3) and we have restricted 
ourselves to consider perturbations of the energy density 
that correspond only to small wavelengths; that is per- 
turbations that can be described by 



p = Po exp 



-2A;|c2|l/2r(l+3a;)/[3(l-Fw)] 



l + 3u} 



Po exp 



1 -I- 3a; 



(4.6) 



Obtaining an exact analytical solution to the differential 
equation (4.7) is very difficult and we will simply consider 
in what follows approximate asymptotic solutions when 
T is very large or very small, for particular values of 
(jj. On the regime of w > —1 and T large, we take the 



observationally plausible value lu = —13/15 so that the 
differential equation (4.7) can then be approximated to 



R"- 



1287rG 
153 



Mo + 



87rGT2 



i?~0. 



(4.7) 



A solution is then again given in terms of the modified 
Bessel function p [32] 



USnGpo, 
135 ' 



(4.6 



For small values of T on the phantom regime at the par- 
ticular case that oj = —5/3, the approximate differential 
equation (4.7) will become 



, 647rG 
R - ^ I MO 



87rGT2 



0, 



(4.9) 



which again admits an approximate solution in terms of 
modified Bessel functions p, that is [32] 



R ~ ri/2 



iMnGpor 



27 



(4.10) 



Choosing now for oq a sufficiently large value then T 
will always be large enough for uj > —1 and small enough 
for LO < —1. Since the most natural and sufficient initial 
condition in the present case is that R « a initially, 
then the only solutions which satisfy the above initial 
condition are for uj = —13/15 (and generally for any uj > 
— 1) the one with p = K, i.e. 



R oc exp I — 




1287rGAto, 



and for u = 
one with p : 



135 



—5/3 (and generally for any u> < 
I, i.e. 



R(xT 



(4.11) 



-1) the 



(4.12) 



Now these two solutions tend to vanish as one lets 
a{t) — > 00, so indicating a complete accelerated gravi- 
tational collapse of the overdense patches, which can be 
modulated by the residual dissipative forces. Although 
the possibility that the black holes that eventually formed 
as the final state of such collapses actually be those su- 
permassive black holes ultimately driving formation of 
the observed galaxies [23,35] is just a conjecture we want 
to advance here (see Sec. V), it already appears inter- 
esting to consider that possibility as a quite promising 
way through which the studied fluctuations may show 
up observationally. 



B. Loss of dark energy 

On the other hand, dark energy does not necessarily 
follow the collapse of dark matter inside the overdensity 
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regions and therefore there could be an energy loss of 
dark energy F in these regions. In the present case this 
energy loss can be described by the equation [36] 



1 




(4.13) 



with R the scale radius of an average overdense region, 

and (j) and V = V{(f)) the scalar field and the potential 
on that region. There are important technical difficulties 
which would appear if we try to calculate F using the 
expressions of R discussed in the precedent subsection. 
In order to evaluate on which cases F ^ we tentatively 
shall simply interpret for a moment the collapsing solu- 
tion ttc = t(i+3")/[3(i+'^)] to the Raychaudhuri equation 
like though if it represented the scale radius of an over- 
density patch, i.e. R = ac, and then use the two ansatze 
considered in Sec. IV for the partition of the full en- 
ergy density and pressure of the unified model into dark 
matter and dark energy components. In such a toy de- 
scription we have for the not partitioned theory 



H 



R 
R 



3lu 



2T 



2{l+tu) 
1 -t- 3a; 



V 



3(l + 3^)-''-(. 



il/2 



327rGT2 



(4.14) 



(4.15) 



(4.16) 



Since no partition of energy density and pressure has 
been made yet, the insertion of Eqs. (4.16)-(4.18) into 
Eq. (4.15) consistently leads to the result that there is 
no loss of any energy on the overdense region and F = 0. 
The use of the first ansatz for partitioning dark energy 
and pressure given by Eq. (3.4) leads to lOj, = — 1, = 0, 
H = -T-i and V = 3/(47rGr2), and hence again F = 0, 
as it was also expected for a cosmological constant. How- 
ever, using the partition given by Eq. (3.6), we obtain 
ujy = 1 — 4>~'^, and finally 



3w f2uj'^ + 2,uj- 1 
~ 2T I l-uj 



(4.17) 



which shows a loss of dark energy in the overdense region 
for any uo > (3+\/T7) /4. We then notice that the collapse 
tachyonic overdense patches will mainly affect the dark 
matter component only if the dark energy component 
does not correspond to a positive cosmological constant 
with (jjy = —1. For in that case the collapse would equally 
affect both components. 



V. A CONJECTURE ON GALAXY 
FORMATION 

Having shown that the evolution of density perturba- 
tions in the unified tachyon theory does not lead to a 



catastrophic growth of perturbations with small wave- 
lengths, but to a negativc-prcssure-assisted super accel- 
erated gravitational collapse of the dark matter and en- 
ergy involved in the overdense regions, leaving rapidly 
formed black holes which may have a large mass, one 
should now consider the observational consequences from 
the formation of such holes. At first sight, the results ob- 
tained in the present paper could in fact be disregarded 
under the claim that no giant black holes of the kind 
discussed in it have been hitherto observed. We are go- 
ing now to conjecture nevertheless that it is precisely a 
family of rapidly forming huge black holes with masses 
analogous to those of the black holes resulting from the 
accelerated collapse of overdense regions in tachyon uni- 
fied dark energy discussed in this paper which have been 
discovered in recent years at the center of most - possibly 
all - quasars (QSOs) and current luminous galaxies [23]. 
The rationale supporting such a conjecture runs up along 
the following points. 

1. The idea that bright galaxies may harbor supermas- 
sive black holes (SMBHs) actually dates back to the year 
1969 and was first advanced by Lynden-Bell [37]. Today 
it is a rather common belief that most - if not all - galax- 
ies have a central SMBH in their bulge [23]. It is also 
currently thought that at the time of the formation of 
one galaxy its present central black hole started also to 
be formed [36]. 

2. The way in which SMBHs were formed is still un- 
clear. They might have a primordial origin, be formed 
from population III stars and their resulting merging, 
or otherwise [38]. The physical process involved could 
respectively be either by directly collapsing from exter- 
nal pressure in the first instant of big bang or by slow 
accreation of matter starting from a stellar size, or oth- 
erwise. The idea that SMBHs may have passed through 
an intermediate observable mass stage along a very long 
accreating process is now gaining increasing support [39]. 

3. The epoch of galaxy formation is thought to be very 
broad, extending from ^; ~ 20 to 2; ~ 2. The galaxies 
were formed from perturbations created at z ^ 1400. 
However, galaxies with redshifts > 2 — 3 are not di- 
rectly observable. Since a baryonic gas inside a dark 
matter potential collapses and forms stars, any result- 
ing massive stars rapidly evolving and fecdbacking heavy 
elements into the gas, the production of first heavy el- 
ements can be approximately equated to the epoch of 
bright galaxy formation. Metal lines have been so far 
identified in QSO spectra [40] indicating that, prior to 
z = 2.5 — 3, the heavy metal line strengths are fairly low, 
with the largest strengths appearing within the interval 
z — 1.5 — 2.5. Thus, non luminous galaxies could have 
existed from nearly the primordial epoch, but the forma- 
tion of a sufficiently high density of stars in them making 
these galaxies bright appears to be fairly recent. 

4. The observational discovery that there exists a sim- 
ple relation between the mass of central SMBHs and the 
speed or the sigma of the stars in the surrounding galaxy 
bulges [41], which was anticipated by the Silk-Rees the- 
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ory [42], seems then to indicate that most of such cen- 
tral SMBHs reached -or were created with - their present 
huge maturity mass just at around z ^ 2 — 3, nearly at, 
or shortly before, the epoch when the giant black holes 
discussed in the present paper probably started to be 
collapsed in an accelerated way. The onset of that epoch 
should nevertheless be marked by the farthest supernova 
at z = 1.7 [43], that is to say around 2 ~ 1. Moreover, 
there exist QSOs at z ^ 6 and bright galaxies at z ^ 3 
which must harbor a SMBH at their central bulge. At 
first sight, these features would prevent SMBHs formed 
by the collapse of overdense tachyon dark energy regions 
to be suitable candidates for being the central holes that 
would induce bulge star creation and hence bright galaxy 
formation. Notwithstanding, there are at least two main 
arguments against that conclusion. On the one hand, it 
has been shown [44] that in the presence of a cosmological 
constant the galaxy formation epoch must have started 
much later and lasted much longer, possibly all the way 
to the present. Similarly, taking into account the Fried- 
mann equation for flat geometry, if instead of the cosmo- 
logical constant density Aq = A/{3Hq) of Ref. [44], we 
introduce the dynamical dark-energy density considered 
in this paper, then a theoretical onset for bright galaxy 
formation can be derived by replacing Ao in the formula 
provided in that reference for the quantity d,TiGp/ {2>Hq). 
Using then Eqs. (2.2), (2.3) and (2.10) one can finally 
have 



1 + 2GF 



H^T^^o sinh^ 



3tGF 

~2T 



-1/3 



(5.1) 



where Hq is the current value of the Hubble constant and 
tcF is the galaxy formation time, tap = td + tc, with td 
the time needed for a galaxy to start forming stars once 
the region with local Hubble parameter H' has collapsed, 
and tc is the collapse time which is in the present case 
determined by 



Jo 



rr 2 



-1/2 



dx. 



(5.2) 



where if the mass of the overdense region of scale radius 
r is denoted by m, 



r{Gm) 



-1/3 



u = 41/3 



and 



u + 



rp2 



= 0, 



(5.3) 



-2/3 



(5.4) 



(5.5) 



with Q, = p/{3H'^), p and H respectively being the mean 
matter density and Hubble parameter of the universe. 



We can check that also in this case for suitable initial 
values of the observable density contrast 5i, the galaxy 
formation epoch is shifted toward later times by some 
two z-units if currently 8ttGp/{3H§) ~ A/(3if^) ~ 0.8, 
while it lasts quite longer, possibly all the way to the 
present, too, such as it occurred in the presence of a 
positive cosmological constant [44] , except in that in the 
present theory zqf depend on time t through the param- 
eter T. It is worth noticing e.g. that in the phantom 
regime where lu < —1, as one approaches the big rip at 
T = 0, — > and zgf — * 0, so that, consistent with 
the ripping off of all materials in the universe, observers 
would detect no galaxies at all. 

On the other hand, and even more importantly, if a di- 
rect interaction term between the two components py and 
Pdm transferring energy from one another is assumed to 
exist and have the form 
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3£p\/QvUjpDM, 



Pv 
3//2' 



SttG 



(5.6) 



then it can be readily seen that the theory discussed in 
this paper becomes completely equivalent to the model of 
coupled dark matter and energy developed by Amendola 
[45]. Therefore one should expect that for our model the 
accelerating expansion would start at a maximum redshif 
of around z ^ 5 [45] in a way which is consistent with 
the type la supernovae Hubble diagram, including the 
farthest known supernova at z = 1.7, while allowing for 
structure formation during the acceleration regime. In 
fact, if there was an epoch at Zacc in the past before which 
uncoupled baryons with critical density ilb dominated a 
decelerating flat universe filled also with tachyon fields, 
then in the present theory with 7 given by Eq. (5.6) that 
epoch would be given by 



1 + 



-(l + 3a;) 



1-^^t 
Oft 



1 l/(3w) 



instead of the currently assumed value [45] 



-(l + 3a;) 



DM 



l/{3u>) 



(5.7) 



1. (5.8) 



Finally, according to Eq. (5.7) for the allowed values of 
Oft, on the phantom regime {ui < —1) the larger \ui\ the 
smaller Zacc- 

5. Thus, since the nature and properties of nonrotating, 
chargeless black holes of any sizes do not qualitatively 
distinguish the original characteristics of the stuff out of 
which they were formed, the difficulties encountered so 
far to justify presently considered mechanisms describing 
how central SMBHs were formed within galaxy bulges 
[46] might all be circumvented if we conjectured that, 
since it turns out that the tachyon field perturbations 
in density contrast took place with a spatial distribution 
nearly matching that of matter fluctuations, and they ini- 
tially may rapidly have grown during a more or less long 
period, shortly after the onset of tachyon dark-energy 
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domination at a sufficiently large z, giant black holes 
mainly made up of dark matter and energy, formed by 
the kind of quickly accelerated collapsing procedure con- 
sidered in Sec. IV, started forming at the galactic centers 
and very rapidly reached a state where they were able to 
switch on stars in the surrounding bulge by means of a 
process similar to the Silk-Rees mechanism [42]. That 
conjecture would seem to decide the debate on which 
came first, the black hole or the galaxy, in favor of the lat- 
ter, but required the existence of galaxies devoid of high 
star densities until their central SMBHs were formed. 

As to the infiuence that ordinary matter sources may 
have on the stability of tachyon dark energy and its 
accelerated collapse, we note that if these sources con- 
tribute as a perturbation, one can analyze the whole 
model via phase space [47]. Then the potentials (2.10) 
and (2.16) and the solutions (2.8) and (2.13) would corre- 
spond to dynamical attractors for, respectively, tachyon 
dark and phantom energy domination with their respec- 
tive cosmological effects. For these attractors we have 
[47] e = W'liy'f = 3 /2, at critical points defined 
by SttGAc which equals 3\/l + w/T for dark energy and 
— 3-\/|cl^| — 1/T for phantom energy. It follows then that: 
(i) the models of Sees. II-IV are stable to the presence 
of small proportions of ordinary matter, and (ii) even 
though that ordinary matter would roughly follow distri- 
bution of the dark matter and energy, one would expect 
it not to significantly undergo initial accelerated collapse 
as it had to obey a different equation of state. Thus, or- 
dinary baryonic matter would essentially remain as the 
stuff out of which stars are going to be formed. 

VI. CONCLUSIONS AND FURTHER 
COMMENTS 

The assumed dark energy content of the universe pos- 
sibly refiects the greatest unsolved problem of all of 
physics. This paper has explored some potentially im- 
portant observational consequences from one of the most 
promising models suggested to describe dark energy, the 
so-called tachyon model. Thus, by using tachyon-like the- 
ories whose starting Lagrangians generally are inspired 
by string theories, we have investigated the properties of 
general FRW cosmological solutions that are fully con- 
sistent with the whole dynamical structure of the theo- 
ries. Assuming a general equation of state p = cop, such 
solutions describe accelerating universes in the interval 
— l/3>a'>— 1 and, among other weird properties char- 
acterizing the phantom regime, all show a big rip singu- 
larity in finite time when lu < —I. It has been generally 
shown as well that the entropy of a universe filled with 
phantom energy is always positive, though its tempera- 
ture becomes definite negative. This keeps still phantom 
cosmology as a real, observationally not excluded possi- 
bility. 

Fluctuations taking place in both of such regimes have 
been studied by using a synchronous perturbation for- 



malism and the relativistic analog of the Newtonian first 
order perturbation theory. In both formalisms the result 
has been obtained that the considered tachyon theory is 
stable to fiuctuations along the entire considered inter- 
val of the state equation parameter —5/3 < u> = p/p < 
-1/3. 

A potential problem with the tachyon models consid- 
ered in this paper stems from the imaginary value of the 
sound speed. In fact a value < implies occurrence of 
instability on scales below the Jeans limit for scalar field 
fluctuations [20] which may grow therefore exponentially. 
Whereas when these models are considered as pure cos- 
mic vacuum components this could in fact be regarded 
as an actual problem, if the tachyon models are viewed 
as the sum of two components [17], one describing the 
negative pressure vacuum stuff and the other describing 
the dust-like cold dark matter contributing ilm ~ 0.35 
and clustering gravitationally at small scales, the results 
obtained in this paper for a negative value for might 
instead be regarded as originating an accelerated pro- 
cess of gravitational collapse inexorably leading to the 
rapid formation of black holes in a way that may ex- 
plain some recent observations in galaxies and superclus- 
ters that concern gravitationally collapsed objects such 
as supermassive black holes. Alcaniz and Lima suggested 
nevertheless that a quintessence component pushes back 
to higher redshift the epoch of formation of some galax- 
ies [48]. The analysis carried out by these authors is 
based however on a simple dark energy component which 
is uncoupled to dark matter. If one assumes that dark 
matter and energy are coupled to each other, then since 
it turns out that the perturbations in density contrast 
may initially grow during a period, we have in fact ad- 
vanced the conjecture that supermassive black holes with 
masses ranging from 10^ to 10^ solar masses which have 
been observed at the center of most galactic bulges could 
actually be originated by the rapid accelerated collapse of 
the gravitationally enlarged overdense regions made up 
of tachyon dark matter and energy. Rather than being 
destructive influences on their surroundings, the result- 
ing giant black holes would then have a very creative 
impact in the formation of the host galaxies. This would 
offer a novel alternate resolution to the problem of su- 
permassive black hole formation. Much calculation and 
observational scrutiny is however needed before such a 
conjecture can be seriously considered. In particular, it 
appears that future observations on supernovae of type 
la at high redshift will be decisive in order to check ac- 
ceptability of our conjecture. 
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